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CONSTRAINTS ON Ua, AND erg FROM GALAXY CLUSTER REDSHIFT DISTRIBUTIONS 
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ABSTRACT 

We show that the number and redshifts distribution of galaxy clusters in future deep cluster surveys can 
place strong constraints on the matter density, Qm, the vacuum energy density, Q\, and the normalization 
of the matter power spectrum, erg. Degeneracies between these parameters are different from those 
in studies of either high-redshift type la Supernovae (SNe), or cosmic microwave background (CMB) 
anisotropics. Using a mass threshold for cluster detection expected to be typical for upcoming Sunyaev- 
Zel'dovich effect (SZE) surveys, we find that constraints on flm and as at the level of roughly 5% or 
better can be expected, assuming redshift information is known at least to z 0.5 and in the absence 
of significant systematic errors. Without information past this redshift, J^a is constrained to 25%. With 
complete redshift information, deep (Mum ~ 10^'^h~^MQ), relatively small solid angle ('^ 12 deg^) 
surveys, can further constrain J^a to an accuracy of ~ 15%, while large solid angle surveys with ground- 
based, large-format bolometer arrays could measure SIa to a precision of ~ 4% or better. 

Subject headings: cosmic microwave background — cosmology: theory — large-scale structure of 
universe — cosmological parameters 



1. INTRODUCTION 

The abundance of galaxy clusters can provide strong 
constraints on cosmological parameters (e.g., Bahcall and 
Fan 1998; Viana and Liddle 1999). Future cluster sur- 
veys, especially those using the Sunyaev-Zeldovich effect 
(SZE), will provide large catalogs of clusters with a se- 
lection function that is remarkably simple (Barbosa et al. 
1996; Holder et al. 2000; Kneissl et al. 2001). Clus- 
ter surveys probe the amplitude of the power spectrum 
as a function of redshift and the cosmic volume per unit 
redshift and solid angle (Haiman, Mohr and Holder 2000, 
hereafter HMH) . This allows unique constraints on cosmol- 
ogy, and provides a sensitive test of structure formation. 
As we show below, the cluster redshift distribution (i.e., 
counts) out to z ~ 0.5 is a powerful probe of and as, 
while a deeper cluster inventory will allow a determination 
of the vacuum energy density, JIa, as well. 

In HMH we showed expected constraints on the equation 
of state of the dark energy and the matter density, assum- 
ing a flat universe and two specific proposed cluster survey 
instruments. In this work, we specialize to the case of vac- 
uum energy (equation of state w = dp /dp = —1), but gen- 
eralize to non-flat universes; we treat the normalization of 
the matter power spectrum erg as a free parameter. Here 
we use less realistic mass limits for cluster detection, to 
focus on the effects of abundance evolution and volume as 
opposed to the particulars of the instruments. Interpreting 
a future cluster survey will require careful consideration of 
the survey selection function. 

In § 2, we outline our methods for calculating survey 
yields, while § 3 outlines our statistical methods. In § 4, 
we present our results and conclude with a discussion of 
possible systematic errors and future work that will be 

^ Hubble Fellow 



necessary to achieve the expected constraints. 

2. CALCULATING CLUSTER SURVEY YIELDS 

To calculate the expected number of clusters per square 
degree as a function of redshift one must consider several 
elements: the comoving volume per unit redshift and unit 
solid angle, dV/dzdfl, the minimum observable mass as a 
function of redshift and cosmology, Miii„(z), and the co- 
moving number density of halos above the mass threshold 
as a function of redshift and cosmology, n{> Mii,„, z). 

The comoving volume per unit redshift is straightfor- 
ward to calculate (e.g., Peebles 1993), and for this work 
we assume that the mass threshold is constant with red- 
shift and cosmology, approximately consistent with ex- 
pectations for surveys using the Sunyaev-Zeldovich effect 
(Holder et al. 2000; Bartlett 2000; Kneissl et al. 2001). 
The mass threshold of detectability in a survey is very 
important, and it is likely that such a mass threshold will 
have a dependence on cosmology. We ignore the cosmolog- 
ical dependence in this work in order to isolate the effect 
of the growth rate of structure as a function of cosmology 
as the primary discriminant between cosmological models. 

For the comoving number density of clusters, we use the 
mass function of dark matter halos obtained from large 
cosmological simulations (Jenkins et al. 2001), where the 
comoving number density of clusters between mass M and 
M + dM is given by 
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where a{M, z) is the variance in the density field smoothed 
on a scale corresponding to mass M and po is the present- 
day mean density of the universe. 

From linear theory, cr(M, z) can be separated into 
a{M,z) = a{M,z = 0)D{z), where D{z) is the linear 
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growth factor. We calculate D{z) directly from linear the- 
ory (e.g., Peebles 1993). The variance a{M,z = 0) is cal- 
culated with a top hat filter of radius appropriate for mass 
M , using the matter power spectrum calculated from the 
fitting functions of Eisenstein and Hu (1999). 

As a fiducial model, we take a flat, low-density model, 
with n,n — 0.3, f^A = 0.7, h — 0.65, erg = 1.0, n = 1, and 
f2{,/i^ = 0.02. We keep h, n and n^h^ fixed for all models 
and vary the remaining three parameters. To isolate the 
effects of the growth of the power spectrum, we artificially 
constrain the power spectrum to be that of the fiducial 
model for all models. In the absence of modifications to the 
power spectrum, these results are absolutely unchanged 
by varying h or flbh^. Relaxing the assumption of a fixed 
power spectrum will slightly affect the constraints on pa- 
rameters, but the qualitative results are unchanged, as the 
shape of the power spectrum has only a mild effect on the 
cluster abundance (HMH) . Large scale structure measure- 
ments, from the spatial distribution of either galaxies or 
galaxy clusters should be able to provide accurate mea- 
surements of the matter power spectrum over the range of 
interest. 

3. ESTIMATING CONSTRAINTS ON COSMOLOGY 

To estimate constraints on cosmological parameters, we 
assume Poisson errors on finely binned (Az = 0.01) model 
predictions. Still finer binning was found to have no effect 
on the results. We determine confidence regions in two 
ways, using a Monte Carlo method and the Fisher matrix 
(Eisenstein, Hu, and Tegmark 1999). 

In the first method, we construct a likelihood space by 
a Monte Carlo method, generating 3000 realizations of the 
input fiducial model, in the form of mock catalogs of clus- 
ters at different redshifts, and then fitting each realization, 
allowing rimj^A and erg to find their best-fit values. In 
this way, we can map out regions of parameter space that 
contain 68% and 95% of the realizations. 

For the case of Gaussian errors, the likelihood is simply 
related to the usual statistic as — 21n£ = x^- In the 
case of Poisson errors, it can be shown that the analogous 
statistic is the Cash C statistic (Cash 1979): 

N 

-2 ln£ = -2^ 71a In Ca - Cq - In na! , (2) 

a=l 

where and are the observed and expected number of 
counts in bin a, respectively, and N is the total number of 
bins. 

A computationally more efficient method to estimate 
confidence regions is to use the Fisher matrix. We con- 
structed the second derivatives of the log likelihood with 
respect to the parameters of interest at the position in 
likelihood space of the fiducial model. Ensemble averag- 
ing then leads to a simple expression for the Fisher matrix: 

-d'^lnC _ deg dcg 1 
dpidpj ^ dpi dpj ma 

where ttIq is the number expected in the fiducial model 
in bin a and the pi are the parameters of interest, in this 
case f2,„, Qjy, and erg. This is also what one would obtain 
assuming a Gaussian probability distribution with a'^ = N 
and using the usual estimator for the likelihood. 



Assuming that the likelihood distribution is approxi- 
mately Gaussian near the peak likelihood, we can use con- 
fidence limits for Gaussian statistics (i.e., x^) to obtain 
68% and 95% confidence regions. 

The advantage of the Fisher matrix method is that it is 
very fast and gives accurate results if the likelihood space 
is well-behaved. While we have no reason to suspect this 
is not the case, we perform the Monte Carlo analysis as a 
useful check. The Fisher matrix approximates confidence 
regions as Gaussian ellipsoids, while the true confidence 
regions could have broad tails or significant curvature. 




Fig. 1. — Top left panel: differential (solid curve) and cu- 
mulative (dashed curve) number counts of clusters as a function 
of redsliift. A mass limit of 10^*/i~^Mq is assumed. Top right 
panel: 68% and 95% confidence levels on Hm-f^A marginalized over 
erg. Shaded regions are from a Monte Carlo analysis, while con- 
tours show results using the Fisher matrix. Dashed contours show 
the current constraints from the Supernova Cosmology Project (see 
http://supernova.lbl.gov), while the diagonal solid line indicates a 
flat universe favored by CMB anisotropics. Bottom right panel; 
Confidence levels on r2m-o"8 marginalized over Qm, similar to the 
top right panel. The diagonal line shows the degeneracy from the 
abundance of massive (Af > 2 X IO^^Mq) local clusters. Bottom 
left panel: Confidence levels on HA-crg (bottom right), marginalized 
over Qm- 

4. RESULTS 

We first look at a relatively low value for the constant 
mass limit, Mn^ = 10^*/i~^Mq, appropriate for a deep 
SZE survey that would cover approximately 12 deg . In 
Figure 1, we show the expected number of clusters as a 
function of redshift, as well as the likelihood contours we 
find in three different projections in the ^m, ^A, cg plane. 
The constraints from the Monte Carlo method are similar 
to those obtained from the Fisher matrix method, espe- 
cially in the Qm-^A plane. This indicates that the likeli- 
hood space is well-behaved, with degeneracies manifesting 
themselves as ellipses. As the Fisher matrix method is 
much faster, this allows a much more extensive study of 
parameters. However, the differences are not negligible for 
the plots involving erg , suggesting that the likelihood space 
in the erg direction is significantly non-Gaussian. 
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Constraints on and ^Ia from cluster surveys are com- 
plementary to other probes. Figure 1 also shows the cur- 
rent constraints on parameters from studies of distant su- 
pernovae, as well as the f2m+f^A=l line expected from in- 
flation and preferred by CMB studies (dc Bcrnardis et al. 
2000; Hanany et al. 2000; Netterfield et al. 2001; Pryke 
et al. 2001). The different orientations of the parame- 
ter degeneracies suggest that joint constraints from CMB 
measurements, supernovae, and deep cluster surveys will 
be a powerful probe of both cosmological parameters and 
potential systematic effects in other data sets. In addition, 
Figure 1 shows that the degeneracy in the Oto-^ s plane is 
in a different direction when compared with local deter- 
minations from massive (M ^ 2 x lO^^h"-'^ Mq) clusters 
(Viana and Liddle 1999). This is due to the different 
mass and redshift ranges probed, with the current local 
determinations coming from much more massive clusters 
(HMH). 

When marginalized over the other two parameters, we 
can expect uncertainties in Q^n, ^a, and as of 6%, 15%, 
and 3%, respectively, assuming complete redshift infor- 
mation. The relatively small number of (~ 400) clusters 
spread over only a few square degrees should allow follow 
up observations to get redshift information. 
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Fig. 2. — Confidence levels on parameters, marginalized over the 
third parameter in the likelihood space, as in Figure 1, but for a 
higher mass limit of 2.5 X 10^^h~^M.(?). Contours show the expected 
68% and 95% confidence regions, using the Fisher matrix. Relative 
axis scaling and the labeling of the curves is the same as Figure 1, 
except in the upper right panel, where we have omitted showing the 
constraints from SNe. Squares, triangles, and X marks show results 
of Monte Carlo investigations of the effects of a systematic tilt in 
the mass function (x cr"'^), a systematic reduction in the mass func- 
tion amplitude (—10%) and a systematic ofi'set in the limiting mass 
(-1-5%), respectively. 

Large ground-based telescopes equipped with bolometer 
arrays will be efficient at detecting clusters through their 
SZE signatures. We assume a typical mass limit for such 
instruments to be 2.5 x lO^'^/f^M© and a typical sky cov- 
erage of 4000deg^ We again calculate constraints on cos- 
mological parameters using the Fisher matrix, with results 
shown in Figure 2. In this case, the possible constraints 
on cosmology are approximately an order of magnitude 



stronger than the first case. The direction of the degener- 
acy is slightly different from Figure 1, due to the different 
redshift distribution and mass range. This indicates the 
potential power of using the distribution of observed clus- 
ters in mass, which we have not utilized in this paper, but 
will explore in future work. The offset between the local 
Oto-cs determination and our central value is an artifact of 
our assumed fiducial model being slightly offset from the 
best fit of Viana and Liddle (1999), although our fiducial 
model is well within their quoted uncertainties. 

Single parameter uncertainties on f^m, Q,a, and erg are 
expected to be approximately 0.7%, 2%, and 0.5% respec- 
tively, for this case, assuming complete redshift informa- 
tion. Again, redshift information will be very important 
for such tests. In this case, redshifts for clusters spread 
over 4000 deg^ is a non-trivial task, but it is clear from 
Figure 2 that the reward for such an exercise will be pre- 
cise measurements of cosmological parameters. 

Also shown in Figure 2 are the effects on parameter es- 
timation of possible systematic effects. The Monte Carlo 
method was used to find the best-fit cosmology when the 
fiducial model was subject to systematic errors. The ef- 
fects of uncertainty in the true mass function were treated 
as either a small tilt, where the mass function has been 
multiplied by cr(M)°-^, or as an overall normalization er- 
ror, where we reduced the amplitude of the mass function 
by 10%. Such effects are allowed within the uncertain- 
ties of Jenkins et al. (2001). The final systematic error 
that we tested was an uncertainty of 5% in the determi- 
nation of the limiting mass. This is significantly smaller 
than current uncertainties in mass estimates, but could be 
achievable with high signal-to-noise X-ray or SZE obser- 
vations of galaxy clusters. Systematic errors could affect 
determinations of erg, with biases in and Oa roughly 
comparable to the la statistical uncertainty. 

We finally adopt a relatively high value for the limit- 
ing mass, A'/iim = 8 x IQ^'^H^^Mq, appropriate for an SZE 
survey like Planck, covering a large fraction of the sky (we 
assume 50%). From Figure 3, it can be seen that such a 
survey would allow strong constraints on cosmological pa- 
rameters, comparable to those expected from the second 
case above. Even with only ~ 1 cluster per square degree, 
the large solid angle allows a large catalog to be assembled. 

In this case, redshift determination would be challeng- 
ing. A reasonable strategy would be to cross-correlate with 
a large solid angle survey such as the SDSS. For example, it 
is reasonable to expect that every cluster found by Planck 
with 2; < 0.5 will appear in the SDSS catalog in regions 
of overlapping coverage. Using only those clusters with 
z < 0.5, and again assuming 50% of the sky is covered 
(although SDSS will only cover ~ 25%), Figure 3 shows 
the expected constraints on parameters. The constraints 
on ag and flm are still strong, but with only the nearby 
sample it is difficult to constrain the vacuum energy. At 
lower redshifts, increasing results in an increased vol- 
ume that largely offsets the effect of the growth function. 
At least in the case that we are considering where the mass 
limits are largely decoupled from cosmology, the clusters 
past z ^ 0.5 are most useful in constraining the cosmologi- 
cal constant. This is unfortunate, as these are the clusters 
that will be most difficult to follow up. 

Single parameter uncertainties on Q^n, and as, as- 
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suming complete redshift information are 1%, 4%, and 
0.9%, respectively for a survey like Planck. For redshift 
information only out to z = 0.5, these numbers become 
2%, 26%, and 2%, respectively. The changes in fi^ and 
as uncertainties are barely larger than would be expected 
from the increased Poisson error from the smaller number 
of clusters. For a ground-based bolometer array, the same 
effect can be observed, with expected uncertainties of 3%, 
23%, and 3%, respectively, without redshift information 
past z — 0.5. 
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Fig. 3. — Confidence levels on parameters, using the Fisher ma- 
trix as in Figure 2, but for a mass limit of 8 X lO^^h-'^MQ. Dashed 
ellipses show the expected 95% confidence regions using only clus- 
ters with z < 0.5, dotted curves show the 95% confidence regions 
using only clusters with z > 0.5, while solid curves show the 95% 
confidence region for the entire sample. Relative axis scaling is the 
same as Figure 1. 

5. CONCLUSIONS AND FUTURE WORK 

We have shown that cluster surveys have the potential 
to be powerful probes of cosmological parameters, and in a 
manner that will serve as a useful check of other methods. 
Efficient follow-up to obtain redshift information will be 
crucial, with most information on JIa coming from clusters 
with z > 0.5. This argument would favor a deep survey on 
a relatively small patch of sky, as follow-up observations 
could be done on a smaller scale. 

The results obtained here are based on highly idealized 
toy models, and are far from realistic in several respects. 
Most importantly, we have not addressed here the effects of 
systematic errors from cluster evolution, which can mimic 
different values of cosmological parameters (but see HMH 
for a quantitative discussion) . In addition, although we ex- 
pect the matter power spectrum to cause relatively little 
cosmological sensitivity (from HMH) , this deserves further 
study We assume that Poisson errors are setting the un- 
certainty on cosmological parameters. With surveys in the 
design stage, it is difhcult to estimate the potential system- 
atic errors that will be important in any interpretation of 
real surveys. 



The mass function, in principle, carries a significant 
amount of information on cosmology and could also be 
a source of systematic error in an interpretation of sur- 
vey results. The mass function has been accurately (to an 
overall error of 30 percent) determined in a handful of cos- 
mologies by Jenkins et al. (2001). It will be important to 
repeat similar numerical experiments in a broader range of 
cosmologies, to have a better understanding of the scaling 
of the mass function with cosmology. Taking full advan- 
tage of the statistical power in the surveys considered here 
will require a description of the mass function, which is 
better than 10% accurate. 

Mass limits as a function of cosmology must be well- 
understood to interpret the source counts. This includes 
both evolution with redshift and scatter in the "detectabil- 
ity" at a given mass. Although temperatures can be ob- 
tained for clusters in a cosmologically independent way 
from their X-ray spectra, the relation of this tempera- 
ture to the halo mass in different cosmologies needs further 
study. Extracting unbiased cosmological parameters from 
real survey data will require multi- wavelength follow-up of 
at least a sub-sample of the detected clusters in order to 
understand the mass limit of the survey to better than 5%. 

Just as large scale structure surveys have learned about 
galaxy evolution and CMB experiments have learned 
about foregrounds, cluster surveys will learn about galaxy 
cluster structure and evolution. As the largest virialized 
objects in the universe, clusters at high redshift will yield 
a wealth of information on galaxy formation and other as- 
pects of cosmology. A study of cluster scaling relations 
that combines SZE and X-ray data can help in separat- 
ing the systematic effects arising from the assumption of 
virialization from cosmology Verde et al. (2001). 

None of these issues should, in principle, prevent robust 
and highly precise constraints on cosmological parameters 
from cluster surveys. The justified use of the Fisher matrix 
will enable fast exploration of many cosmological parame- 
ters. The small subset explored here and in HMH suggest 
that these constraints should be highly complementary to 
large scale structure studies (e.g., Phillips et al. 2000), 
CMB anisotropies, and probes of distance, such as type la 
supernovae. In addition, tracing the redshift evolution of 
the cluster abundance offers a unique opportunity to track 
the growth of structure that is expected from hierarchical 
clustering, providing a strong test of the underlying cold 
dark matter model. 
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